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ABSTRACT 

The bright Type Il-plateau supernova (SN) 2004dj occurred within the young, massive stellar cluster 
Sandage-96 in a spiral arm of NGC 2403. New multi- wavelength observations obtained with several 
ground-based and space-based telescopes are combined to study the radiation from Sandage-96 af- 
ter SN 2004dj faded away. Sandage-96 started to dominate the flux in the optical bands starting 
September 2006 (~800 d after explosion). The optical fluxes are equal to the pre-explosion ones 
within the observational uncertainties. An optical Keck spectrum obtained ~900 d after explosion 
shows the dominant blue continuum from the cluster stars shortward of 6000 A as well as strong 
SN nebular emission lines redward. The integrated spectral energy distribution (SED) of the cluster 
has been extended into the ultraviolet region by archival XMM-Newton and new Swift observations, 
and compared with theoretical models. The outer parts of the cluster have been resolved by the 
Hubble Space Telescope, allowing the construction of a color-magnitude diagram. The fitting of the 
cluster SED with theoretical isochrones results in cluster ages distributed between 10 and 40 Myr, 
depending on the assumed metallicity and the theoretical model family. The isochrone fitting of the 
color-magnitude diagrams indicates that the resolved part of the cluster consists of stars having a 
bimodal age distribution: a younger population at ^10-16 Myr, and an older one at ~32-100 Myr. 
The older population has an age distribution similar to that of the other nearby field stars. This 
may be explained with the hypothesis that the outskirts of Sandage-96 are contaminated by stars 
captured from the field during cluster formation. The young age of Sandage-96 and the comparison 
of its pre- and post-explosion SEDs suggest 12 < A/ prog < 20 M Q as the most probable mass range 
for the progenitor of SN 2004dj. This is consistent with, but perhaps slightly higher than, most of the 
other Type Il-plateau SN progenitor masses determined so far. 

Subject headings: supernovae: individual (SN 2004dj) — galaxies: individual (NGC 2403) 
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1. INTRODUCTION 

The theory of stellar evolution predicts that massive 
stars (M > 8 Mq) end their lives as core-collapse su- 
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pernovae (CC SNe, e.g.. IWooslev et al.| [2002h In par- 
ticular, after the main-sequence phase the most massive 
stars undergo heavy mass loss, become stripped stellar 
cor es and explode a s Type Ib/c supernovae (SNe Ib/c; 
see lFilippenko|[l997l for a discussion of SN spectral clas- 
sification). Stars close to the lower mass limit of CC are 
thought to produce Type Il-plateau SNe (SNe II-P). Re- 
cent observations to detect the progenitors of CC SNe 
support this scenario. Currently there are 10 SNe II 
(1987A, 1993J, 1999ev, 2003gd, 2004A, 2004et, 2005cs, 
2006ov, 2008bk; the last seven are Type II-P) whose pro- 
genitors have been directly identified in pre-explosion im- 
ages ( seelHendrv et al1l2006HLi et alJ2007HMattila et all 
l2008tlLeonard et al.ll2008HSmartt et al.ll2008l and refer- 
ences therein), and the mass estimates are M < 15- 
20 M Q for all of them. Moreover, upper mass lim- 
its were derived for a number o f other SNe II from 
nondetections of t heir p r ogenitors (IVan Dvk et al.1 120031 : 
iMaund fc Smarttl l2005t iLeonard et alj I2008D . and the 
highest upper limit was found to be M 20 M©. These 
observations have led to the conclusion that SNe II- 
P likely originate fro m "low-mass" massive progenitors 
with M < 20 M Q (|Li et al.l I2006L 120071 : iSmartt et all 
l2008h . and the fate of stars with M > 20 M Q may be 
a SN Ib/c explosion. 

On the other hand, the progenitors of SNe Ib/c (even 
the brightes t and closest ones) hav e escaped direct detec- 
tion so far (|Crockett et al.|[2007af ). The most promising 
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candidate is SN 2007gr, which occurred in a compact, 
massive stellar cluster in NGC 1058 that has been de- 
tected wit h the Hubble Space Telescope (HST) prior to 
explosion (jCrockett et alj|2008h . 

SN 2004dj, the closest (~ 3.5 Mpc. IVinko "eTaT]l2"00ll 
and one of the brightest SNe since SN 1987A, was a 
Type II-P event. It occurred within a young, massive 
cluster, Sandage-96 (S96) in NGC 2403. SN 2004dj has 
been extensivel y studied th r ough multi- wavelength ob- 
servations; see I Vinko et al.l (|2006( ) (hereafter Paper I) 
for references. In particular, several attempts were 
made to infer the mass of the progenitor by compar- 
ing the pre-explosion magnitudes and colors of S96 
with theoretical spectral energy distributions (SEDs) 
to determine the age, hence the turnoff mass, of the 
cluster. These resulted in a range of possible pro- 
genitor masses from M « 12 M to M > 20 M Q 
dependi ng on the assumed meta ll icity and/o r red- 
dening dMaiz-Apellaniz et alj 120041: I Wang et al.l 120051 : 
IVinko et al.l 120061 ) . However, all of these studies suf- 
fered f rom the age-reddening an d age-metallicity degen- 
eracy (R enzini &: Buzzonil 119861 ) , because the available 
pre-explosion observations covered only the optical and 
near-infrared (NIR) bands. 

The main aim of this paper is to derive further con- 
straints on the progenitor mass of SN 2004dj from post- 
explosion observations of S96, made after its reappear- 
ance from the dimming light of the SN. The cluster 
has been successfully redetected both by our ground- 
based and space-based observations, showing no signifi- 
cant change in its optical light level with respect to the 
pre-explosion level. We have extended the wavelength 
coverage of the observed SED into the ultraviolet (UV) 
with new Swift/ WOT observations. Moreover, the clus- 
ter has been partially resolved by our new HST/ ACS 
observations, which provides an additional opportunity 
to infer age constraints on its stellar population. 

The paper is organized as follows. In §2 we present 
our new multi-wavelength observations made from the 
ground and space. The SED of the cluster and its color- 
magnitude diagram are constructed and fitted with the- 
oretical model predictions in §3. We discuss the results 
in §4 and present our conclusions in §5. 

2. OBSERVATIONS 

2.1. Optical Data 

2.1.1. Ground-Based Photometry 

The light variation of SN 2004dj in the nebular phase 
was followed from Konkoly Observatory (see Paper I for a 
description of the telescopes and detectors). In addition, 
Johnson B, Johnson V, and narrow-band Ha images 
were taken with the 90Prime camera on the 2.3 m Bok 
telescope at the Steward Observatory, Arizona (Fig. [l] 
left panel). 

The magnitudes of SN 2004dj were calculated via aper- 
ture photometry based on the same sequence of local 
standard stars as in Paper I. The photometric data ob- 
tained after 2005 May are summarized in Table [TJ 

The light curves are plotted in Figure [51 Following 
the usual decline in the nebular phase (starting ~ 100 d 
after explosion), the light curves approached a constant 
level around day 800; see the two right-hand panels of 
Figure [U 



As expected, the flattening of the light curves is caused 
by the increasing contribution of the radiation from S96, 
emerging from the fading light of SN 2004dj. In Fig- 
ure [2] the dotted horizontal lines mark the pre-explosion 
magnitudes of S96 (see Paper I). 

From the two right-hand panels of Figure [2] it is appar- 
ent that the post-explosion magnitudes of S96 are almost 
identical to the pre-explosion ones in V , R, and /. There 
is a very slight excess in the B band (~ 0.1 mag), which 
is about the same as the photometric uncertainty of the 
data. Although it cannot be ruled out that this excess is 
due to some kind of systematic error in the calibration of 
the -B-band data (the deviation from the pre-explosion 
level is ~ la), it is interesting that the ground-based B 
and V magnitudes agree very well with those obtained 
by Swift/WOT (see §2.2.2). 

2.1.2. Keck Spectroscopy 

A late-time spectrum of SN 2004dj (exposure time 
2200 s) was obtained on 2006 Dec. 23 (~ 900 d after 
explo sion) with the DEIMOS spectrograph (iFaber et al.l 
2003) mounted on the 10 m Keck-II telescope in Hawaii. 
The 1200 line mm -1 grating was used, with a slit 
1.1" wide, resulting in a resolution (full width at half- 
maximum intensity) of 2.7 A. The slit w as aligned close 
to the parallactic angle (Filippenko 1982), so differential 
light losses were not a problem. 

As seen in Figure El the Keck spectrum is clearly a 
composite of S96 and the nebular ejecta of SN 2004dj. 
Longward of 6000 A, strong emission lines of Ha, [O I] 
AA 6300, 6363, and [Fe II] A7155 A, characteristic of a 
typical nebular SN II-P spectrum at late phases, can be 
identified. Shortward of 6000 A, the blue continuum 
dominates the spectrum; Na I D appears in emission, 
which emerges mostly from the SN ejecta, but H/3 is in 
absorption. Clearly, the radiation from the young stellar 
population of S96 is visible in this regime. The shape of 
the spectrum is fully consistent with the predictions of 
population-synthesis models (see Paper I and §3). 

2.1.3. HST Observations 

SN 2004dj and its surrounding area was observed with 
HST/ ACS on 2005 August 28, -425 d after explosion 
(GO-10607; PI. B. Sugerman). Four sets of 4 drizzled 
frames were obtained through the F606W and F81AW 
filters, and 3 sets were recorded with the ^43511^ filter. 
In the latter case the POLUV polarization filter was also 
placed in the beam. This made it possible to study the 
polarization of the SN light, but slightly complicated the 
photometry of the F435W frames, causing a systematic 
shift of the zcropoint in the standard transformation (see 
§3.2). 

The ACS frames, reduced and calibrated by the HST 
pipeline (including MultiDrizzle), were downloaded from 
the HST archive at the Canadian Astronomy Data Cen- 
tre 16 . Because SN 2004dj was still very bright compared 
with the rest of S96 at the epoch of these observations, 
its point-spread function (PSF) was subtracted from the 
ACS frames. We used the TinyTim software 17 (version 
6.3) for calculating the ACS PSFs in each filter. Since the 

16 CADC is operated by the National Research Council of 
Canada with the support of the Canadian Space Agency. 

17 http://www.stsci.edu/software/tinytim/tinytim.html . 
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Fig. 1. — Left panel: color-combined B, V, and Ho frames of NGC 2403 obtained with the 2.3 m Bok telescope at Steward Observatory 
on 2007 Jan. 28 (UT dates are used throughout this paper). Right panel: Swift/XJVOT image of NGC 2403 obtained on 2007 Dec. 2 (it, 
uvwl, and uvw2 filters were selected as red, green, and blue channels). The field of view of both frames is 10.0' X 7.2'; north is up and east 
is to the left. The position of S96/SN 2004dj is marked. 



TABLE 1 

Late-Time BVRI Photometry of SN 2004DJ. 



UT Date 



JD - 
2,450,000 



^expl 

(d) 



B 

(mag) 



17.79 (0.08) 
17.89 (0.11) 
17.93 (0.11) 
17.98 (0.07) 
18.22 (0.10) 
18.01 (0.09) 
18.15 (0.09) 
18.11 (0.10) 
18.15 (0.06) 
18.13 (0.08) 
18.11 (0.08) 



V 
(mag) 



17.25 (0.03) 
17.68 (0.05) 
17.71 (0.05) 

17.85 (0.03) 
17.79 (0.04) 
17.83 (0.04) 
17.88 (0.04) 

17.86 (0.04) 

17.86 (0.02) 
17.76 (0.04) 

17.87 (0.03) 



R 
(mag) 



16.69 (0.08) 

17.19 (0.11) 

17.38 (0.11) 

17.50 (0.07) 
17.53 (0.10) 

17.53 (0.09) 

17.54 (0.09) 

17.51 (0.10) 

17.51 (0.06) 
17.47 (0.06) 



I 

(mag) 



Instrument 



16.44 (0.06) 
16.86 (0.09) 

17.00 (0.09) 

17.07 (0.06) 

17.08 (0.08) 
16.96 (0.08) 
17.08 (0.08) 
17.08 (0.08) 

17.01 (0.08) 
16.99 (0.08) 



2005- 
2006- 
2006- 
2006- 
2006- 
2006- 
2006- 
2006- 
2007- 
2007- 
2007- 



11- 09 
01-27 

08- 23 

09- 07 

09- 22 

10- 17 

12- 22 
12-27 

01- 28 

02- 09 

03- 06 



3684.6 
3762.5 
3971.6 
3986.3 
4001.6 
4026.6 
4092.5 
4097.6 
4128.0 
4141.4 
4166.3 



500 
577 
787 
801 
817 
842 
907 
913 
943 
956 
981 



Konkoly 0. 

Konkoly 0. 

Konkoly 0. 
Konkoly 

Konkoly 0. 

Konkoly 0. 

Konkoly 0. 

Konkoly 0. 
Steward 
Konkoly 
Konkoly 



6 m Schmidt 
6 m Schmidt 
6 m Schmidt 
1.0 m RCC 
6 m Schmidt 
6 m Schmidt 
6 m Schmidt 
6 m Schmidt 
2.3 m Bok 
1.0 m RCC 
1.0 m RCC 
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Fig. 2. — BVRI light curves of SN 2004dj from ground-based photometry. The horizontal lines mark the pre-explosion magnitudes of 
S96. In the left panel the scaling on the abscissa is logarithmic. The two right-hand panels show the same data as the left-hand ones, but 
focus on the region around 800 d. 

analytical PSF works less effectively for drizzled frames, 
the flatfield-corrected ".FLT" frames were used for the 
PSF removal. After sub-pixel registration, the individ- 
ual frames belonging to the same filter were averaged. 
The model PSF was then scaled to the peak of the SN 
and subtracted from the combined frame. 

The result is shown in Figure |U The encircled region 
(r = 35 pixels ss 15 pc) contains S96 with its unresolved 
inner and resolved outer parts. Several bright red and 
blue giants are visible in the outer region. The color of 



the unresolved inner part is also very blue, in accord with 
the pre-explosion photometric observations and the pro- 
posed young cluster age (see Paper I) . It is also apparent 
that SN 2004dj occurred near the projected center of the 
cluster (there are some artifacts due to the incomplete 
PSF removal at the SN position, but they are less than 
1% of the subtracted SN flux). 

Photometry of the stars appearing in the ACS 
frames was obtained with the DOLPHOT 18 1.0 software 

18 http://purcell.as.arizona.edu/dolphot . 
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4500 5000 5500 6000 6500 7000 7500 
Observed wavelength (A) 

Fig. 3.— Optical spectrum of SN 2004dj/S96 obtained with the 
10 m Keck-II telescope on 2006 Dec. 23. The identified bright 
emission lines are formed in the SN ejecta. 




Fig. 4. — Color-combined image from the HST F435W (blue), 
F606W (green), and F814W (red) frames of S96 obtained on 2005 
Aug. 28. The field of view is about 4" X 4"; north is up, east to 
the left. The green circle has a radius of 35 pixels (~ 15 pc at the 
distance of NGC 2403), and it defines the boundary of the cluster 
region in our analysis. The PSF of SN 2004dj was modeled and 
removed (see text). 

(jDolphinl 120001 ). DOLPHOT incorporates corrections 
for geometric distortions of the ACS camera, cosmic- 
ray removal, object identification, PSF fitting (using pre- 
computed PSFs via TinyTim), charge-transfer efficiency 
correction, and transformation into standard photomet- 
ric systems. It works best with the flatfield-corrected 
".FLT" frames. All of these frames were processed with 
DOLPHOT, and the resulting magnitudes belonging to 
the same filter were combined frame by frame. Only 
those stars that could be identified on at least two frames 
with the same filter were retained in the final list. The 
photometric errors were computed from the scatter of 
the individual data around their mean value, taking 
into account the individual magnitude errors computed 
by DOLPHOT. The final magnitudes were converted to 
Johnson- Cous i ns B, V, and / using the calibration by 
iSirianni et al.l (|2005( ). Note that 0.3 mag has been added 
to the transformed B magnitudes to take into account 
the transmission of the HST POLUV polarization filter, 
which was used together with the filter during 

the observations. The results are analyzed in §3.2. 

2.2. Ultraviolet Data 



2.2.1. XMM-Newton Observations 

Prior to the explosion of SN 2004dj , S96 was observed 
with the Optic al/UV Monitor tele scope (OM) on board 
XMM-Newton (jMason et al.ll200lD on 2003 April 30 (PI. 
M. Pakull). The FITS frames and tables containing the 
photometric data (reduced and calibrated by the SAS 
pipeline) were downloaded from the XMM-Newton Sci- 
ence Archive 19 . The instrumental magnitudes of S96 (ob- 
ject #1057) are listed in Table [5] Unfortunately, no B 
or V observations were made, so full transformation into 
the standard Johnson system cannot be computed. How- 
ever, applying the UV transformation equations in the 
OM Calibration Documentation 20 , the correction in the 
U band is only 0.019 mag; thus, the instrumental mag- 
nitudes in Table [2] should well represent the Vega-based 
standard magnitudes of S96. 

Finally, the observed count rates were transformed into 
fluxes (in erg s _1 cm~ 2 A -1 ) using the conversion fac- 
tors listed in the OM Calibration Documentation (CAL- 
TN-0019-3-2, p. 17, Table 4). It is known that such a 
conversion is only approximate, because it depends on 
the SED of the object. However, the SED of S96 in the 
blue/UV regime is very similar to that of an early- type 
star or a white dwarf (cf. §3.1). The count rate to flux 
conversion was calibrated using white dwarf spectropho- 
tometric standards. Therefore, the flux conversion of S96 
is fairly robust. 

2.2.2. Swift Observations 

The Swift Observatory (jGehrels et all |2004[ ) was 
launched into orbit on 2004 N ovember 20. Its Ultr avio- 
let/Optical Telescope (UVOT. Homing et al.l [20051 ) was 
used to observe SN 2004dj/S96 at five epochs. Table [3] 
summarizes the basic parameters of these observations. 
A color-combined UV image (made from the data ob- 
tained on 2007 December 3) is presented in the right- 
hand panel of Figure [TJ 

The UVOT observations were downloaded from the 
Swift data archive 21 . The sky-subtracted frames were 
processed and analyzed with a self-developed script in 
the following way. First, the individual exposures be- 
longing to the same filter (stored as extensions of the 
same FITS file) were co-added with the UVOTIMSUM 
routine of the HEAsoft software 22 . For those frames that 
were obtained without binning, the summed frames were 
rcbinned by 2 x 2 binning (1" x 1") in order to enhance 
the signal-to-noise ratio (S/N) of the point sources. The 
fluxes of S96 and the local photometric standard stars 
(that could be identified on the UVOT frames) were com- 
puted with aperture photometry in IRAF 23 . The photo- 
metric calibr ation was do n e acco rding to the latest pre- 
scriptions by iPoole et al.1 |2008). The aperture radius 
was set as 5" (5 pixels on the 2x2 binned frames), while 
the sky was computed as the "mode" of the pixel values 

19 http://xmm.esac.esa.int/xsa/index.shtml . 

20 http:/ /xmm. vilspa.esa.es/external/xmm_sw_cal/calib/index.shtml 

21 http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl . 

22 http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ . 

23 IRAF, the Image Reduction and Analysis Facility, is dis- 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 



Sandage-96 After the Explosion of SN 2004dj 



5 



TABLE 2 

XMM-Newton OM Observations of S96 



Obs. ID 


UT Date 


Exp. time 


UVW2 


Flux" 


UVWl 


Flux" 


u 


Flux" 






M 


(mag) 




(mag) 




(mag) 




0150651101 


2003-04-30 


6304 


16.87 


9.20 


16.76 


6.59 


17.37 


4.04 






a 


0.11 


0.93 


0.06 


0.37 


0.06 


0.23 



Note. — "The flux units are 10 16 erg s 1 cm 2 A 1 . 



TABLE 3 

Swift UVOT Observations of SN 2004dj/S96 



Obs. ID 


UT Date 


Exp. time 


uvw2 


Flux" 


uvm2 


Flux' 1 


uvwl 


Flux" 


u 


Flux' 1 


b 


Flux" 


V 


Flux" 






w 


(mag) 




(mag) 




(mag) 




(mag) 




(mag) 




(mag) 




00035870002 


2006-10-09 


2215 










16.96 


6.98 


17.44 


3.45 


18.02 


3.59 


17.65 


3.24 


00035870003 


2006-10-15 


5553 










16.91 


7.33 


17.25 


4.08 


17.94 


3.88 


17.81 


2.81 


00035870004 


2007-04-06 


2308 


16.84 


9.61 


16.91 


6.92 


17.00 


6.78 


17.13 


4.58 










00036563001 


2007-12-03 


6440 


16.81 


9.89 


16.80 


7.66 


16.95 


7.09 


17.31 


3.88 


18.24 


2.93 


17.76 


2.94 


00036563002 


2007-12-06 


2382 


16.99 


8.32 


17.01 


6.31 


17.07 


6.30 


17.24 


4.13 


18.07 


3.44 


17.80 


2.83 






average 


16.88 


9.27 


16.91 


6.96 


16.98 


6.90 


17.27 


4.02 


18.07 


3.46 


17.76 


2.96 






a 


0.10 


0.84 


0.11 


0.68 


0.06 


0.39 


0.11 


0.41 


0.13 


0.40 


0.07 


0.20 



Note. — "The flux units are 10 16 erg s 1 cm 2 A 1 . 



in an annulus with r„ = 10 and r out = 15 pixels cen- 
tered on the point source. The summed, sky-corrected 
fluxes in analog-to-digital units (ADU) were divided by 
the dead-time corrected exposure time (defined by the 
keyword EXPOSURE in the FITS headers) to obtain 
the raw count rates in ADU s _1 . These raw count rates 
were c orrected for coincidence loss following iPoole et al.l 

Hop). 

Finally, the corrected count rates were transformed 
into magnitudes and physical fluxes by using t he pho- 
tometric calibrations given bv IPoole et all (|2008H . (Note 
that we have applied the formulae based on the Pickles 
stellar spectra, instead of gamma-ray burst spectra, be- 
cause the SED of S96 is more like that of a star than 
a gamma-ray burst.) No color-term correction was ap- 
plied to the magnitudes, since we intend to compare 
physical fluxes rather than magnitudes from UVOT and 
other instruments. The color-term corrections have been 
computed only for checking the deviation of the UVOT 
magnitudes from the magni tudes in the stan dard John- 
son/Bessell system for S96 (IPoole et al.l [2008). The cor- 
rections are U — u = 0.22 mag, B — b = 0.03 mag, and 
V — v = 0.03 mag, where lower-case letters refer to the 
Swift filters. It is seen that the UVOT b and v magni- 
tudes are fairly close to the standard system, while the u 
magnitudes are slightly brighter. The final UVOT fluxes 
and their uncertainties are listed in Table [3j they are 
analyzed further in §3.1. 

From our ground-based optical photometry (§2.1) it 
was concluded that SN 2004dj faded below the light level 
of S96 by September 2006 (~800 d after explosion). Be- 
cause the UV flux of SNe II-P diminishes more rapidly 
than the optical flux (jlmmler et alj l2007t iBrown et al.l 
I2007D . the UVOT observations made in October 2006 
and April 2007 recorded mostly the cluster light. Fig- 
ure [5] shows the comparison of the UV fluxes with the 
Keck spectrum (dotted line) and ground-based fluxes (as- 
terisks). The UV fluxes observed by different satellites 
before and after explosion agree nicely within the errors. 
The agreement is also very good in the B and V bands, 
between the ground-based and space-based observations; 




-16 ' — 1 ' 1 1 ' 1 1 1 

2000 3000 4000 5000 6000 7000 8000 9000 
Wavelength (A) 

Fig. 5. — Observed SED of S96 in the UV and optical region 
plotted against the observed wavelength. The observatories are 
indicated in the legend. 

the differences are ~ 0.06 mag and ~ 0.02 mag in B and 
V, respectively. Similar agreement was obtained for some 
of our local photometric standard stars, although most 
of them were too bright for UVOT photometry. 

3. RESULTS 

We now present the analysis of the observations de- 
scribed in the previous section. 

3.1. SED Fitting 

The physical properties of S96 were discu ssed by 
iMaiz-Apellaniz et all (|2004h . iWang et alj (|2005h . and in 
Paper I. These studies were based on pre-explosion pho- 
tometry of S96 (broad-band Johnson UBVRI and the 
14-color BATC system in the optical, and 2MASS JHK 
in the NIR). Fitting the optical through NIR SED with 
single stellar population (SSP) models, all of these stud- 
ies revealed that S96 is a young, compact stellar cluster 
with an age of ~8-20 Myr. The uncertainty is caused 
by the strong age-reddening and age-metallicity correla- 
tions in the SED fitting, and also to the sensitivity of 
the stellar-evolution models applied for constructing the 
SED of an SSP with a given age. 

Extending the wavelength coverage of the observed 
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TABLE 4 

Time-Averaged SED of S96 Before and After SN 2004dj. 
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Fig. 6. — Comparison of pre- and post-explosion SEDs of S96. 

SED may help break the age-reddening-metallicity de - 
generacy (R cnzim fc Buzzoni|[l986l ; iKavirai et al.ll2007l ). 
The SED of young stellar clusters can be best charac- 
terized in the UV, because the UV luminosity, origi- 
nating from the most massive, fast-evol ving supergiants , 
strongly correlates w ith the cluster age (|0'Connelllll999t 
iBuzzoni et alJ l2Q07f ) . By adding the UV data from 
XMM-Newton and Swift (§2.2) to the optical through 
NIR SED used in the previous studies, one can get a 
better constraint for the cluster age, and hence the SN 
mass. 

There is, however, an additional, non-negligible source 
of systematic uncertainty in the interpretation of the 
SEDs of young massive clusters: they contain ~ 10 4 - 
10 5 stars, and the high-mass end of their initial mass 
function (IMF) is poorly populated. Because these 
stars are also the most luminous ones, the observed 
SEDs of clusters having the same physical parameters 
will statistically deviate from the model SED depending 
on the actual number of t h eir most massive, luminous 
stars (ICervino et all l2002b ICervino k. Luridianal I2004 
Damet et all l2004t ICervino fc L uridiana 2006|). This 
clearly limits the applicability of model SEDs computed 
by assuming an analytical, completely sampled IMF (i.e., 
practically an infinite number of stars). Because this ef- 
fect was not taken into account in the previous studies 
of S96, here we address it in detail. 

In Table H] we have collected all available photometric 
data (converted to fluxes in erg s _1 cm~ 2 A -1 ) of S96, 
including both pre-explosion and post-explosion observa- 
tions. We have used our BVRI photometry made 800 d 
after explosion (Table [T]) for the representation of the 
post-explosion flux in the optical. Unfortunately, there 
are no post-explosion observations in the JHK bands at 
our disposal. 

The pre- and post-explosion SEDs are plotted together 
in Figure [S] It is apparent that the two datasets agree 
within the errors. The agreement in the UV and optical 
bands implies that the removal of the flux of the progen- 
itor of SN 2004dj from the integrated light of the cluster 
has not caused a significant loss of light in these bands. 

In order to fit theoretical SEDs to the observations, 
we have defined an averaged "normal" SED of S96 by 
combining the pre- and post-explosion data. In the UV 
range, between 2000 and 4000 A, we adopted the average 
of the fluxes from XMM-Newton/ OM and Swift/VYOT. 
The ground-based data in this spectral range are ex- 
pected to be less reliable than the satellite-based ones, 
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2025 
2825 
3360 
3450 
3663 
4210 
4361 
4540 
4872 
4925 
5270 
5448 
5795 
6075 
6282 
6407 
6656 
7057 
7546 
7776 
7980 
8023 
8480 
9182 
9739 
12200 
16300 
21900 



450 
750 
360 
900 
650 
320 
890 
340 
1280 
390 
340 
840 
310 
310 
1150 
1580 
480 
300 
330 
1230 
1540 
260 
180 
260 
270 
2130 
3070 
3900 



9.20 
6.59 
3.11 
4.04 
3.58 
3.28 
3.33 
2.88 
2.74 
2.54 
2.22 
2.72 
2.18 
2.14 
1.93 
1.75 
1.90 
1.98 
1.79 
1.64 
1.37 
1.59 
1.50 
1.35 
1.24 
1.06 
0.70 
0.28 



After • 



uvw2 
uvm2 
uvwl 



B 
V 
R 
I 



2030 
2231 
2634 
3501 
4329 
5402 
4361 
5448 
6407 
7980 



760 
510 
700 
875 
980 
750 
890 
840 
1580 
1540 



9.27 
6.96 
6.90 
4.02 
3.46 
2.96 
3.75 
2.69 
1.78 
1.39 



0.93 
0.37 
0.62 
0.23 
0.43 
0.12 
0.29 
0.08 
0.10 
0.16 
0.10 
0.11 
0.08 
0.11 
0.10 
0.15 
0.10 
0.16 
0.09 
0.10 
0.17 
0.11 
0.23 
0.13 
0.25 
0.08 
0.07 
0.05 



plosion 



0.84 
0.68 
0.39 
0.41 
0.40 
0.20 
0.34 
0.25 
0.23 
0.24 



This paper (XMM) 
This paper (XMM) 
Wang et al. (20051 
This paper (XMM) 
Larsen (1999) 
Wang et al. (20051 
Paper I 

Wang et al. (20051 
Davidge (2007^ " 
Wang et al~T2005l 

Wang et al. (20051 
Paper I 

Wang et al. (20051 
Wang et al. (20051 
Davidge (20071 
Paper I 

Wang et al. (20051 
Wang et al. (20051 
Wang et al. (20051 
Davidge (20071 
Paper I 

Wang et al. (20051 
Wang et al. 120051 
Wang et al. 120051 
Wang et al. (2005) 
Skrutskie et al. (1997) 
Skrutskie et al. (19971 
Skrutskie et al. ( 1997) 



This 
This 
This 
This 
This 
This 
This 
This 
This 
This 



paper 
paper 
paper 
paper 
paper 
paper 
paper 
paper 
paper 
paper 



(Swift) 
(Swift) 
(Swift) 
(Swift) 
(Swift) 
(Swift) 



Note. — A c and AA denote the central wavelength and the full 
width at half-maximum intensity (FWHM) of a given filter, in A. 
The flux units are 10 -16 erg s _1 cm -2 A -1 . Uncertainties are 
given in parentheses. 



because of the higher probability of systematic errors in- 
troduced by the local atmospheric conditions. In the 
optical we used the Johnson-Cousins BVRI data. In 
the NIR, only the p re-explosion 2MASS JHK fluxes 
(|Skrutskie et al.lfl997t ) were available to us. 

There is a possibility that the observed SED is some- 
what contaminated by foreground/background stars be- 
longing to NGC 2403, altering the fluxes from being en- 
tirely due to a SSP. The amount of this contamination is 
difficult to estimate, because S96 itself may contain some 
older stars captured from its galactic neighborhood (see 
§3.2 and §4). However, the background subtraction we 
applied during photometry should have removed most of 
the flux from background stars. Due to the compact- 
ness of S96 (< 30 pc diameter, Figf?]), the number of 
foreground field stars should be minimal (see §3.2 for 
a detailed discussion) . Consequently, the SED fluxes are 
expected to be due mostly to S96. Since S96 is by far the 
brightest source in this region, the contamination from 
field stars should not exceed the estimated errors of the 
"normal" SED fluxes, ~ 10%. 

During the analysis the distance of NGC 2403 was fixed 
at 3.5 Mpc. This optimal value is found by combining 
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various distance measurement results for SN 2004dj and 
its host galaxy, as discussed in Paper I. 

In order to test whether the effect of statistical 
IMF sampling allows the m odeling of the cluster SED, 
ICervino fc Luridianal j2004) introduced the concept of 
the "lowest luminosity limit" (LLL). It can be simply 
expressed as follows: the integrated luminosity of the 
cluster must be higher than that of the most luminous 
star of the model isochrone, at any wavelengths. The 
LLL is a strong function of age and it also depends on 
the considered wavelength regime (or filter band) . As ex- 
pected, LLL gives the strongest constraint for the m odels 
with age 1-100 Myr (jCerviho fc Luridianal 12004). Be- 
cause the possible ages of S96 found in previous studies 
are all in this interval, it is necessary to check whether 
the cluster meets this criterion. The normal SED fluxes 
described above were dereddened and converted to ab- 
solute magnitudes assuming a distance of 3.5 Mpc and 
E(B - V) = 0.07 mag (see Paper I). This resulted in 
Mb — 9.2 mag, Mr ks —10.0 mag and increasing up to 
Mk ~ —12.3 mag. The predictions for the most massi ve 
stars from the Padova isochrones (jCioni et al.l l2006aH bT) 
with t = 8 Myr (the youngest age proposed for S96 so far) 
are -7.3, -9.1, and -11.0 mag (Z = 0.008) and -8.2, 
-9.1, and -10.1 mag (Z = 0.019) for the B, R, and 
K bands, respectively. We see that S96 is at least ~ 1 
mag brighter in the optical/NIR bands, and it is also ~ 2 
mag brighter in U . This criterion becomes more relaxed 
toward higher ages as the possible most luminous stars 
become fainter. Our first conclusion is that although S96 
is close to the LLL for 8 Myr, it is definitely above it, 
so this cluster is sufficiently rich to make the comparison 
of its observed SED with theoretical models statistically 
feasible. 

For comparison with the observations, we have applied 
three different classes of SSP models with somewhat dif- 
ferent input physics. First, as in Paper I, we sel ected the 
GALAXEV models of lBruzual fc Chariot! (|2003D that are 
based on the Padova evoluti onary tracks. Secon d, we 
chose the SPEED models by Uimenez et~aH (|2004D that 
were computed using a new set of stellar interior mod- 
els, evolutionary tracks, and different treatment of the 
mass loss. For both of these models, a Salpeter IMF 
Was adopted, similarly to Paper I and previous stud- 
ies. The machine-readable data of these two model sets 
were downloaded from the SPEED website 24 . Third, we 
appli ed the SSP models generat ed by the Starburst99 25 
code ([Vazquez fc Leithererll2005l ). The Starburst99 mod- 
els are highly configurable; the user may choose between 
different evolutionary tracks, atmospheric models, and 
pre-computed spectral libraries to create a unique set 
of SSP models. In order to test the model dependency 
of the results, we have chosen the Geneva evolutionary 
tracks and Kroupa IMF, and generated SSP SEDs using 
metallicities Z = 0.004, 0.008, and 0.02, between t = 
and 100 Myr. Note that the metallicity resolution of the 
SPEED models is lower; only models with Z = 0.004 
and 0.02 are available. The age step of the Starburst99 
models was selected as At = 1 Myr sampling linearly 
between 1 and 100 Myr, providing much better age res- 

24 http:/ /www. astro. princeton.edu/~raulj/SPEED/index. html 

25 http://www.stsci.edu/science/starburst99/ . 



TABLE 5 

Parameters of the SED Fitting, Without IMF 
Fluctuations (See text) 



Model 


Z 


T c 


M c 


E(B - V) 


x 1 






(10 6 yr) 


(10 3 M Q ) 


(mag) 




J04 


0.004 


8 


27 


0.09 


1.422 


J04 


0.020 


21 


99 


0.04 


2.313 


BC03 


0.004 


35 


114 


0.13 


2.794 


BC03 


0.020 


26 


92 


0.13 


0.536 


BC03 


0.020 


9 


37 


0.17 


0.886 


SB99 


0.004 


10 


90 


0.08 


3.279 


SB99 


0.008 


9 


26 


0.12 


3.521 


SB99 


0.020 


9 


24 


0.10 


1.624 


SB99 


0.020 


40 


91 


0.07 


3.918 



olution than the SPEED models for t > 10 Myr. 

The model SEDs were compared with the observations 
via the usual x 2 fitting. Note that throughout the paper 
we use the reduced x 2 (the sum of the squares of residuals 
divided by the number of data points). The optimized 
parameters were the cluster mass M c , the cluster age 
T c , and the reddening E(B — V), while the metallicity 
and the distance were kept fixed. The interstellar extinc- 
tion at any wavelength in the considered UVOIR spec- 
tral regime was calculated adoptin g Ry = 3.1 and the 
averag e Galactic reddening law of iFitzpatrick fc Massal 
(12001 . 

First, the fitting was computed in the canonical way, 
not taking into account the statistical fluctuations in the 
IMF sampling, and the models were fitted directly to the 
observed normal SED, as if the cluster were composed of 
an infinite number of stars. In this case the \ 2 function 
was defined as 

1 N oba 

X 2 = ^^4^° 6 W-^c-S mod (Ai,r c )] 2 , (1) 

Nobs ^ *i 

where N b s is the number of observed points in the nor- 
mal SED, F obs is the observed flux at wavelength (cor- 
rected for distance and extinction), Oi is its uncertainty, 
and S mod {\i,T c ) is the flux of the model SED with age 
T c , at the same wavelength. Because the model fluxes 
are usually normalized to 1 M , the cluster mass M c en- 
ters simply as a scale parameter in this expression. The 
results of these calculations are summarized in Table [5j 
The age resolution for the SPEED models was slightly 
increased by interpolating between the two neighboring 
model SEDs for those ages that were not covered by the 
original models. However, without this correction, the 
best-fitting parameters did not change significantly. The 
cluster parameters inferred in this way are very similar 
to the results of earlier investigations cited above. 

Second, the statistical IMF sampling was taken into ac- 
count as follows. For each age and wavelength, an uncer- 
tainty cr mo d{Xi,T c ) was assigned for any model flux as a 
measure of the fluctuation of the model SED fluxes due to 
the random sampling of the IMF. Then, each model flux 
was modified as S(Aj)±£ where < £ < cr mo d(K,T c ) is a 
random variable. This step was repeated N mo d (= 1000) 
times, thus constructing a series of model SEDs that 
fluctuate around the original model fluxes. For the fcth 
model (1 < k < N mo d), x\ was computed as in Eq. (1). 
Finally, following the recommendation by an anonymous 
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Fig. 7. — The relative SED flux uncertainty due to IMF fluctua- 
tions as a function of age. Different symbols are used for different 
filter bands, as indicated by the labels in the upper-right corner. 

referee, the final \ 2 was determined as 

x 2 = -^(^ISrX ( 2 ) 

where Py. = exp(— 0.5N o f, s xl) is the likelihood that the 
fcth model describes the observations. Eq. (2) means 
that the Pk likelihoods are averaged, and the final x 2 
is computed from P ave . This approach gives lower final 
X 2 values than the simple average of the individual x\ 
values. 

The modification of the x 2 function in Eq. (2) ensures 
that the \ 2 is mostly sensitive to those models that are 
particularly affected by the sampling effect (i.e., those 
whose <J mo d is high) but give a good fit to the obser- 
vations, while giving lower weight to those models that 
produce inferior fits. Also, the x 2 value remains mostly 
unchanged when the random sampling effect is negligi- 
ble, because in this case the individual \\ values (and the 
corresponding -P^s) are nearly the same for each random 
model. 

Of course, the reliability of this approach depends 
heavily on the proper selection of the <J mo d values. More- 
over, the <J mo d values belonging to different filters are cor- 
related, because the addition or subtraction of one bright 
star would affect the cluster flux in all bands. This corre- 
lation is not reflected by our random models, as the fluc- 
tuations were added to the fluxes independently in each 
band. The proper treatment of all these effects would 
require full computation of detailed SSP models, taking 
into account the statistical population of the IMF, which 
is, however, beyond the scope of this paper. 

Instead, we followed a simpler approach as a first ap- 
proximation: a Salpeter IMF between 0.5 and 100 M Q 
was randomly populated with stars until the sum of their 
mass reached M c = 10 5 M Q (roughly the mass of S96). 
The luminosities of these stars in the UBVRIJHK 
bands were selected from Padova isochrones at a given 
age (between 4 and 100 Myr) and simply summed up to 
get an estimate of the cluster SED. One thousand such 
SEDs with the same cluster mass and age were generated 
and the standard deviations of the fluxes were computed 
at each wavelength. 

In Figure [Jj the relative flux uncertainty (defined as 
Cmod/S, where S is the mean flux of the SED at a given 
wavelength and <r mo d is its standard deviation) for each 
band is plotted as a function of age for Z = 0.019. The 



results are very similar in the case of Z = 0.008 and Z = 
0.004. It is apparent that the relative uncertainty drops 
below 5% for t > 10 Myr, but it may reach up to 30% for 
younger clusters at certain bands. Generally, the U and 
B bands are less affected (the maximum uncertainty is 
<10%), but the scatter becomes higher toward the NIR. 
This means that the SED of such young clusters would 
fluctuate from cluster to cluster with amplitudes up to 
30% of the flux values. The reason for the drop of the 
fluctuations above 10 Myr is that the IMF fluctuations 
are most pronounced at M > 20 M Q for a 10 5 M cluster 
considered here, and the highest possible mass decreases 
quickly below this limit for t > 10 Myr. 

The relative flux uncertainties (a m od / S (Xi)) are listed 
in Table [6l These values were applied in the computa- 
tions of the modified x 2 , as discussed above. The results 
of the x 2 minimizations are in Table [JJ Note that in 
this case the interpolations between ages were not ap- 
plied at all, and the fitting was computed only for those 
ages that were covered by the original models. However, 
when we used the simple average of the x\ values as the 
final x 2 ■> instead of the Pk likelihoods as in Eq. (2), the 
parameters of the best-fitting models did not change. 

Figure [5] shows the results of the SED fitting with the 
IMF sampling effects taken into account (note that only 
the mean fluxes of the best-fitting models SEDs are plot- 
ted). The fitting of the unperturbed SEDs (i.e., ignoring 
the IMF fluctuations) resulted in very similar figures. In 
general, the models applied in this study give an ade- 
quate representation of the observed SED of S96 with 
the cluster parameters collected in Tables [S] and [7J 

In many cases, the x 2 map (upper-right panel of Fig. [8]) 
showed not one, but two distinct minima at two different 
ages, regardless of the presence or absence of IMF fluc- 
tuatio ns. This was first noted by i Maiz-A pellaniz ct al. 
(|2004l) , and it is confirmed here. IMaiz-A pellani z et al. 
(|2004j ) found that their younger solution (T c w 14 Myr) 
had the lower x 2 of the two minima. In the present case, 
it turned out to be model dependent. In the case of the 
BC03 models, the older solution has slightly lower x 2 , 
while for the SB99 models it is the younger one that has 
a deeper minimum. 

Figure [HI shows the % 2 of the best-fitting models (i.e., 
those listed in Tables [5] and [7J plotted as a function of the 
cluster age T c . It is apparent that all the models with 
the lowest x 2 have solar metallicity [Z = 0.02) except 
one model by J04 which has Z — 0.004. Also, the canon- 
ical model fits (without IMF fluctuations, open symbols) 
clearly show the two-age structure noted previously: a 
"young" solution with T c ss 8-10 Myr and an "older" 
one with T c « 25-40 Myr, preferring the "young" solu- 
tions on the basis of % 2 - On the other hand, the models 
with IMF fluctuations (filled symbols) present a continu- 
ous distribution between the same age limits, also with a 
preference for the younger, solar- metallicity models. The 
disappearance of the bimodal distribution is mainly due 
to the fact that when random IMF sampling is taken 
into account, the younger, less massive models are more 
affected, which may increase their x 2 ■ (Though there 
might be a few models that fit the observations very well, 
this is probably not true for the majority of them.) This 
gives higher preference to the models that are less af- 
fected by random sampling, thus shifting the best-fitting 
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TABLE 6 
Relative Flux Uncertainties ov, 



,odl &mod{^i) Tc). a 



Xgi 

< 10 Myr 
> 10 Myr 



uvwl 
0.05 
0.03 



uvm2 
0.05 
0.03 



uvwl 
0.05 
0.03 



0.08 
0.03 



0.10 
0.03 



o.io 

0.03 



0.15 
0.03 



— r~ 

0.20 
0.03 



—J— 

0.25 
0.03 



0.30 
0.03 



0.30 
0.30 



Note. — a The values for the UV bands are based on extrapolation, since these bands were not covered by the applied isochroncs. 
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Fig. 8. — SEDs of SSP models fitted to the observations of S96 with IMF sampling effects taken into account. Filled circles represent the 
"normal" SED of the cluster constructed by averaging the observations within the selected wavelength bands (see text). The horizontal 
error bars denote the wavelength range of the broad-band filters, while the vertical error bars correspond to the uncertainties of the data. 
Solid and dotted lines show the mean fluxes of the theoretical SEDs with parameters indicated as labels. The top-right panel presents the 
grayscale map of a \ 2 space of a bimodal solution, where the dark areas represent the most probable solutions (the horizontal stripes are 
due to the finite resolution of the grid). 



TABLE 7 
Parameters of the SED Fitting 



Model 


Z 


T c 


M c 


E(B - V) 


x 1 






(10 6 yr) 


(io 3 Me) 


(mag) 




J04 


0.004 


8 


34 


0.11 


1.802 


J04 


0.02 


20 


94 


0.04 


2.511 


BC03 


0.004 


40 


121 


0.10 


2.973 


BC03 


0.02 


20 


8(5 


0.14 


0.629 


BC03 


0.02 


8 


34 


0.18 


1.453 


BC03 


0.02 


10 


41 


0.15 


1.492 


SB99 


0.004 


30 


72 


0.09 


3.257 


SB99 


0.008 


11 


30 


0.08 


3.018 


SB99 


0.02 


14 


39 


0.11 


1.497 


SB99 


0.02 


28 


72 


0.09 


2.666 


Note. - 


- "IMF fluctuations taken into account (see 


text). 



models toward those with T c > 10 Myr. 

Although the lowest-x 2 models have Z = 0.02, this 
may be misleading, because the metallicity is a weakly 
constrained parameter in SED fitting. Additional in- 
formation about the possible cluster metallicity may 
help strengthen the confidence of this parameter. In- 



Z=0.020 
Z=0.008 
Z=0.004 



▲ 



A 
▲ 



O 

10 



O 



□ 



• o 

20 30 40 

Age [Myr] 



50 



Fig. 9.— The x 2 of the best-fitting model SEDs as a function of 
cluster age (in Myr). Filled symbols denote the models with IMF 
fluctuations (Table |7j, while open ones indicate the canonical solu- 
tions (without random IMF population, Table[5]l. Different symbol 
shapes denote different metallicities as indicated in the upper-right 



dependent observations suggest that the average metal- 
licity of NGC 2403 is below solar. The oxygen abun- 
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dance in NGC 2403 at the position of S 96 is [O/H] = 
—0.2 4, from spectroscopy of H II regions (Pilvugi n et al.l 
120041) . The distribution of red supergiants in the 
color-magnitud e diagram of stars within the inner disk 
(Davidge 2007) also suggests that the average metallic- 
ity of the whole population is Z » 0.008. On the other 
hand, iMaiz-Apellaniz et al.l (|2004l ) preferred solar metal- 
licity for S96 based on its galactocentric dist ance and the 
abund ance gradient in NGC 2403 found bv lFierro et alj 
l(l986l) . There seems to be no clear consensus on the pos- 
sible metallicity of S96. The fitting of the cluster SED 
as a whole suggests Z as 0.02, but with high uncertainty. 
This problem is investigated further in the next section 
with the analysis of the resolved stellar population of 
S96. 

Contrary to previous solutions (Paper I) , the reddening 
parameter is much more tightly constrained in this case. 
Its average value is E(B — V) « 0.10 ± 0.05 mag, while in 
previous studies values as high as E(B — V) ~ 0.35 mag 
were also proposed. Also, the age-reddening-metallicity 
degeneracy is much reduced in this case, owing to the 
increased wavelength coverage of the observed SED in 
the UV. This reddening value is in very good agreement 
with E{B-V)sn = 0.07±0.1 mag derived for SN 2004dj 
(Paper I). Although the SED of the whole cluster can, 
in principle, be much more affected by intracluster red- 
dening than SN 2004dj itself if the position of the SN 
within the cluster is on the near side toward the ob- 
server, our new HST observations also strongly suggest 
that E(B — V) w 0.1 mag for all the stars resolved within 
the cluster (§3.2). We adopt E(B - V) = 0.1 ± 0.05 mag 
for the rest of this paper. 

Looking for additional constraints on the cluster pa- 
rameters, we attempted to fit the high-resolution Keck 
spectr um (§2.1.2) with high-resolutio n SSP model spec- 
tra by Gonzalez Dclgad o~et al.l (|2005h that are based on 
the Geneva tracks. The motivation for this was the sen- 
sitivity of s ome spectral featur es on the cluster age and 
metallicity (|Koleva et al.ll2008t ). However, this analysis 
was complicated by the obvious presence of the SN neb- 
ular lines that dominate the red part of the integrated 
spectrum. Because H/3 is in absorption, its origin should 
be mostly from cluster stars, but the contribution from 
SN 2004dj may be non-negligible. It is difficult to esti- 
mate the SN contamination at ~ 900 d, because very few 
observed SN spectra exist at this epoch. We examined 
three spectra o f SN 1987A. tha t were taken 700-1000 d 
past explosion (|Pun et al.lll995h . and found that there is 
an emission feature that may be attributed to H/3 in these 
spectra. The amplitude ratio of these emission compo- 
nents, Ha/H/3, is found to be ~ 4. 

The SN 2004dj contamination at H/3 in the Keck 
spectrum was estimated in the following way. First, a 
Lorentzian emission profile was fitted to the observed 
Ha line, taking into account the absorption component 
of the other cluster stars from the SSP models (note that 
the Ha emission due to SN 2004dj is so strong that the 
absorption component has only a minor effect on the fit- 
ted amplitude). Second, this profile was shifted to the 
rest wavelength of H/3 with its amplitude divided by the 
Ha/H/3 amplitude ratio (AR, using AR= 4 as default) 
and its damping parameter 7 multiplied by 0.548 (tak- 
ing into account that the Lorentzian FWHM for pressure 
broadening scales with A 2 ). The fitting was recomputed 



using different amplitude ratios between 3 and 6, to test 
the sensitivity of the results on this parameter. 

The fitting was computed so that the sum of the 
SSP model spectrum and the adopted H/3 emission pro- 
file with fixed AR and 7 was fitted to the observed 
spectrum by varying the age and the cluster mass. 
E(B — V) was fixed at 0.1 mag and the model spectra 
were reddened using the s ame Galactic reddening law 
(|Fitzpatrick & Massal [2007h as for the SED fitting. We 
restricted the computation of % 2 in the vicinity of the 
H/3 line, between 4820 A and 4910 A rest wavelengths. 

It was found that the H/3 profile can be fitted satis- 
factorily with a broad range of ages, depending on the 
chosen metallicity and Ha/H/3 AR. Figure fTUl shows two 
of the best-fitting models with Z = 0.02 and AR= 4 
(left figure) and 6 (right figure). The corresponding 
cluster ages are 25 Myr and 9 Myr, respectively. As- 
suming AR= 4, the ages of the best-fitting models were 
found to be between 25 and 40 Myr depending on metal- 
licity. However, they turned out to be 30 - 45 Myr 
for AR= 3, and 9 - 25 Myr for AR= 6. The results 
from SSP models based on Pa dova evolutionary tracks 
([Gonzalez Delgado et al.ll2005h showed very similar be- 
havior, but resulted in higher cluster ages and poorer fits 
(i.e. higher \ 2 )- 

It is concluded that in general, the fitting of the high- 
resolution spectrum confirms that the cluster is probably 
younger than 50 Myr, but the H/3 profile turned out to 
be mostly sensitive to the non-negligible contamination 
from SN 2004dj. As a result, the line-profile analysis 
could not lead to a unique solution for the cluster age 
and metallicity. Hence, at present, we cannot use the 
H/3 line-profile analysis to further constrain the cluster 
age. 

We have also examined the hypothesis that S96 may 
not be a SSP resulting from a single, rapid initial star- 
burst. Although a single starburst is a more plausi- 
ble mechanism for the formation of a massive, compact 
stellar cluster, continuous star formation is t aking place 
within the disk of NGC 2403 (lDavidgdl200l . We have 
checked whether the SED of S96 could be fitted by that 
of a SSP resulting from continuous star-formation rate 
(SFR); Starburst99 models were compared with Geneva 
tracks, Kroupa IMF, and different metallicities assuming 
continuous SFR. The SFR was simply scaled to match 
the U-band observed flux of the cluster SED. Two of the 
models with Z — 0.008 metallicity are plotted in Fig- 
ure HU Regardless of age, these models are too bright in 
the UV and too faint in the NIR, which suggests that 
the observed SED cannot be described by continuous 
SFR. The same result has also been obtained using other 
metallicities, or applying the Padova evolutionary tracks. 
Note that the presence of a hypothetical dense, intraclus- 
ter dust cloud may significantly alter the shape of the re- 
sulting SED, but a detailed study of such a model would 
require much better observational coverage of S96 at IR 
wavelengths. 

3.2. Isochrone Fitting 

The computed photometry of the H ST/ ACS frames 
(§2.1.3) was used to construct color- magnitude diagrams 
(CMDs) of S96 using either B - V or V - I as color. 
We have selected and examined all resolved stars within 
R = 35 pixels (~ 15 pc) around the cluster center (green 
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Fig. 10. — The fitting of the deredshifted Keck spectrum with high-resolution SSP models having Z = 0.02 (sec text). Each panel shows 
the observed Keck spectrum (blue), the best- fitting SSP spectrum (green), and the best-fitting SSP spectrum with the assumed Ha and 
H/3 SN emission lines added (red). The upper, smaller panels zoom in on the H/3 and Ha regions. The assumed amplitude ratio is indicated 
at the top of the figures. 

region are indeed physically associated with S96, and not 
just a random concentration of unrelated field stars. 

The separation of the cluster members and the field 
stars can also be illustrated in their magnitude his- 
togram. In Figure fl3l the relative frequency (i.e., the 
number of stars in a magnitude bin divided by their total 
number) of the field stars (filled bars) and those within 
the cluster area (open bars) that have V — I > 1 mag 
is plotted as a function of the observed V magnitude. 
The distribution of these red stars clearly indicates that 
in the cluster area there is a significant excess of stars at 

V ~ 22. 5 ±0.5 mag. Their magnitude distribution can be 
roughly approximated by a Gaussian, and it is markedly 
different from that of the field stars, being monotonically 
increasing toward fainter magnitudes. 

Turning back to Figu re [T2") it also cont ains the lat- 
est Padova isochrones (|Cioni et al.l 2006a. b) including 
variable molecular opacities in the thermally pulsing 
asymptotic giant branch (TP-AGB) phase, assuming so- 
lar metallicity. The ages of the plotted isochrones (10, 
16, 32, 63, and 100 Myr) are indicated in the leg- 
end. The isochrones were reddened with E(B — V) = 
0.1 mag (§3.1) assuming the Galactic reddening law 
(|Fitzpatrick k. Massd[2007l) . and shifted to the 3.5 Mpc 
distance of the host galaxy. This reddening value seems 
to be a good estimate for the other field stars as well. 
The Z — 0.019 tracks were selected, because the fitting 
of the integrated cluster SED produced the best results 
using this metallicity (see §3.1). Comparing the CMDs 
with isochrones of Z = 0.008 and 0.004, it was found 
that these isochrones do not extend enough to the red (to 

V — I ~ 2 mag) where some of the bright cluster stars 
reside. However, the age distribution of the observed 
stars (i.e., the concentration of stars along the computed 
isochrones) is the same as in the case of Z = 0.019, so 
the age limits of the resolved population of S96 are found 
to be rather insensitive to the actual metallicity of the 
cluster. 

It is interesting that in the CMDs the field stars follow 
roughly the same distribution as the cluster members. 
Note that the blueward distribution of all stars in the V 
vs. B — V diagram below 23 mag, and the redward dis- 
tribution in the other CMD below 23.5 mag, are due to 
the decreasing sensitivity of the detector/filter combina- 
tion in that color regime (that is, the incomplete detec- 



Wavelength (A) 

Fig. 11. — Comparison of continuous SFR models with the ob- 
served SED of S96. The models were generated by the Starburst99 
code based on Padova AGB-enhanced tracks with Z = 0.008 and 
Kroupa IMF. The ages of the models are indicated in the legend. 

circle in Fig. [3| as possible cluster members. Note that 
the visible diameter of the unresolved inner part of the 
cluster is ~ 15 pixels, corresponding to ~ 6 pc at the 
distance of NGC 2403. 

The CMDs are plotted in Figure 021 where the filled 
circles denote the possible cluster members, within the 
R = 35 pixels radius (referred to as the "cluster region" 
hereafter), while crosses represents the other field stars, 
outside the cluster region. 

Thirty stars have a measured V — I color within the 
cluster region. There are 21 such stars with B — V color. 
However, only 7 stars are common to the two samples, 
due to the reduced sensitivity of ACS in the blue. 

The field-star contamination within this region was es- 
timated by putting outside the cluster region an annulus 
having the same area as that of the cluster region, and 
counting the stars within this annulus. Using different 
inner radii for the annulus, but keeping its area fixed, 
the number of field stars was found to vary between 1 
and 5. Adopting its mean value, the expected number of 
field stars within the cluster region is 3 ± 2. The relative 
contamination of projected field stars within the cluster 
region is ~ 10%. Assuming that the positions of field 
stars follow a Poisson distribution with A = 3 as the ex- 
pected value, the probability of the occurrence of 8 field 
stars within the cluster region (i.e., ~ 26% contamina- 
tion) is ~ 0.8%. This number strongly suggests a 99% 
probability that at least 22 stars found within the cluster 
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Fig. 12. — HST color-magnitude diagrams of S96 using B — V (left panel) and V 
cluster stars (within the encircled region in Fig. 5), while the plus signs denote field stars. In the right panel, those cluster stars that have 
both B — V and V — I colors are plotted as triangles. The fitted Padova isochrones for Z = 0.019 are also shown (labels indicate the cluster 
age in Myr). 
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Fig. 13. — Relative frequency of field stars (filled bars) and 
cluster stars (open bars) with V — I > 1 mag as a function of 
observed V magnitude. The dotted curve is a Gaussian illustrating 
the presence of the bright cluster stars at V ~ 22.5 ± 0.5 mag. 

tion of objects). The completeness limit (the magnitude 
limit above which all stars are detected regardless of their 
color) was estimated as V ~ 22.5 (My ps —5.2) mag for 
the V vs. B-V diagram, and V » 23.5 (My « -4.2) 
mag for the V vs. V — I diagram. In order to have better 
statistics, in the following we analyze the V vs. V — I 
diagram. 

From Figure [T2l the age of each cluster star was deter- 
mined as the age of the nearest isochrone. In some cases, 
when different isochrones ran very close to each other, 
only upper and lower limits (e.g., 63 < t < 100 Myr) 
could be determined. 

It is apparent that the brightest cluster stars are closest 
to the ~ 10 Myr isochrone consistently in both diagrams. 
However, there are only 2 or 3 such stars, so they may 
also be binaries consisting of older/fainter stars. Most 
of the bright resolved stars have V — I « 2 mag and 
are distributed between the 10 and 16 Myr isochrones. 
These are in very good agreement with the ages of the 
SED fitting with the lowest x 2 (§3-1). Because these stars 
are expected to have the most significant contribution to 



32 
Age [Myr] 

Fig. 14. — The histogram of the computed ages of cluster stars. 



the integrated cluster SED, this agreement gives further 
credibility to the age estimates found in §3.1. 

On the other hand, 16 cluster stars out of 30 (^50% of 
the resolved cluster population) are close to or below the 
32 Myr isochrone. The detection becomes increasingly 
color dependent below 23 mag, so the actual number of 
such stars may be higher. There are a few very red stars 
at ^25 mag, where we cut the observed sample, because 
the errors calculated by DOLPHOT started to exceed 1 
mag (note that the real brightness uncertainties of these 
stars may be higher, but we used the errors given by 
DOLPHOT as a selection criterion). 

FigHH shows the histogram of the ages (the age res- 
olution follows that of the isochrones). About half of 
the resolved cluster members fall into the 10-16 Myr age 
interval, while the other 50% have ages distributed be- 
tween 16 and 100 Myr. These results suggest that the 
resolved population of S96 cannot be represented by a 
single age. Instead, a "young" population with an age of 
10-16 Myr and an "old" population at 30-100 Myr seem 
to exist within the cluster area. 

It is interesting to compare the spatial distribution of 



Sandage-96 After the Explosion of SN 2004dj 



13 



450 
460 
470 
480 
490 
500 
510 
520 



560 570 580 590 600 610 
X (pixels) 



620 



Fig. 15. — The spatial distribution of "young" stars (filled circles) 
and "old" stars (open symbols) within the cluster region. The 
approximate center of the cluster is marked with a "+" sign. 

the "young" and "old" stars in the cluster area. This 
is shown in Figure 1151 where the image coordinates of 
the resolved stars (in pixels) are plotted. Filled circles 
represent the "young" stars, and open symbols denote 
the "old" ones. The "young" (bright) stars appear to be 
concentrated around the cluster center, while the "old" 
stars are more scattered. This is illustrated further in 
Figure [16j where the V magnitudes and the V — I colors 
are plotted as a function of the radial distance from the 
cluster center (in pixels). Again, the "young" stars seem 
to dominate the inner area within r ~ 6 pixels. On 
the other hand, the light from the unresolved part of the 
cluster is also strong here, making the detection of fainter 
stars very difficult in the central area. Thus, the lack of 
"old" stars in this area is surely affected by selection. The 
right-hand panel of Figure [16j suggests that the average 
color of the "young" population is somewhat redder than 
that of the "old" population. 

Summarizing the results obtained in this section, we 
conclude that the resolved stars in S96 belong to two pop- 
ulations having ages of 10-16 Myr and 30-100 Myr. The 
"young" stars are brighter and somewhat redder that the 
"old" ones, and they are located closer to the central part 
of the cluster. The age limits of the "young" population 
are in good agreement with the lowest % 2 models in the 
SED fitting (§3.1), when the random population of the 
IMF was taken into account, but slightly higher than the 
age of the best- fitting canonical models (~ 8 Myr), i.e., 
those without random IMF population. 

3.3. The Absence of Ha Emission Around S96 

NG C 2403 is know n to show intense star-forming ac- 
tivi ty (iDav idgc 2003). From deep gri and JHK imag- 
ing, [Davidgc (2007) found that the SFR during the past 
10 Myr has been ~ 1 M© yr _1 in the whole disk of 
NGC 2403. The SFR was highest in the region at galac- 
tocentric distances 2-4 kpc. The intense star formation 
in the inner disk may explain the existence of young 
(~ 8-10 Myr) compact clusters, such as S96 which is 
at R G c ~ 2.7 kpc. 

Young clusters are able to ionize the surrounding hy- 
drogen clouds, showing up as large, bright H II re- 
gions. The measured Ha luminosity i s known to cor- 
relate with the SFR of these co mplexes (Kcnnicu 
iPflamm-Altenburg et~aT1 120071 ). The ionizing UV pho- 
tons come mostly from the young, massive OB stars lo- 
cated inside the clouds. Because the lifetime of such 



stars is short, the number of ionizing photons decreases 
rapidly after ~7-8 Myr for clusters /asso ciations that 
were formed after an initial starburst (jDopita et al.l 
2006). Thus, the presence/absence of Ha emission 
around S96 may give an additional, independent con- 
straint on the age of the cluster. 

Figure [T7] shows the color-combined image of 
NGC 2403 obtained with the 2.3 m Bok telescope at 
Steward Observatory (see §2) using B, V, and Ha fil- 
ters for the blue, green, and red colors, respectively. It 
is apparent that there are a number of extended H II 
regions showing Ha emission in the vicinity of S96 (the 
marked object), as expected in a stellar field with ongoing 
star formatio n . Foll owing the method applied recently by 
iRamva et al.l (|2007j ), the SFRs of these complexes were 
estimated to be 0.01-0.001 M Q yr _1 , typical of such Ha- 
emitting regions. However, S96 appears stellar, without 
any indication for extended Ha emission. This suggests 
that the flux at Ha is coming entirely from inside the un- 
resolved cluster. Indeed, it is very likely that the source 
of this emission is mostly from SN 2004dj (§2.1.2). 

The lack of any extended Ha emission around S96 can 
be used to estimate a lower limit for the cluster age, as 
outlined above. The number of ionizing UV photons as 
a function of age was estimated by the Starburst99 code 
(see §3.1) applying Geneva tracks, Salpeter IMF (but 
neglecting random IMF sampling), and Z = 0.02. The 
cluster mass was fixed at M c — 50, 000 M Q , between the 
cluster masses derived during the SED fitting (see Ta- 
ble [5j . The calculated numbers of ionizing photons have 
been converted to the radius of the H II region applying 
the formula 

where Q(H°) is the number of photons capable of ion- 
izing hydrogen, N e is the number density of electrons 
(complete ionization was assumed: N e = N p « Nh), 
and qg is the eff ective recombination coefficient for H 
(|Osterbrock| [l989). The value of ctB was estimated using 



2.5!)K)- |i (^; 



-0.833 



(4) 



assuming T e = 10 4 K (|Moore et al.|[2002T) . 

In Figure [TS] the radius of the ionization zone is plot- 
ted as a function of the cluster age. The continuous line 
shows the results for N e = 100 cm -3 (a typical electron 
density in bright H II regions), while the dashed and 
dotted lines illustrate the results if N e was an order of 
magnitude higher or lower. Note that changing the clus- 
ter metallicity down to Z — 0.004 caused only negligible 
alterations in these curves. It is apparent that at ~ 10 
Myr the ionized cloud has ~ 5 pc radius, which is simi- 
lar to the radius of S96 as seen by HST/ACS (§2.1.3). 
Above 10 Myr the radius quickly decreases. At ~ 20 Myr 
it is only ~ 1 pc, which is much less than the size of the 
cluster. This suggests that the 10 < t c < 20 Myr cluster 
age found in the previous sections is consistent with the 
lack of resolved H II region around S96. 

For an age of ~ 8 Myr, which was proposed by the fit- 
ting of SEDs without random IMF fluctuations, n ~ 
10 pc, which is slightly larger than the visible cluster size. 
Thus, the ~ 8 Myr age may be less probable than the 
t c > 10 Myr ages found above. However, if N e > 100 
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Fig. 19. — Time dependence of masses of stars with M > 7 Mq 
from Padova tracks. Labels indicate the mass in Mq. Each curve 
ends at the age of the last model of the corresponding track, indi- 
cating the lifetime of a star with the given initial mass. 



Fig. 17.— False-color image of NGC 2403 in the vicinity of S96 
(red, Ha; green, V; blue, B) obtained with the 2.3 m Bok telescope 
at Steward Observatory. The field of view is about 2' X 2' ; north is 
up and east to the left. The marked object is S96. The extended 
red areas are nearby H II regions. 
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Fig. 18. — The calculated radius of the ionization zone (in pc) as 
a function of cluster age (in Myr). The continuous line corresponds 
to N e = 100 cm -3 , while the dashed and dotted lines illustrate the 
dependence of the result on this parameter. 

cm~ 3 is allowed, i?H n can be easily reduced to ~ 5 pc 
at - 8 Myr. 

It is concluded that using the absence of an Ha- 
emitting region around S96 results in a lower limit of 
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Fig. 20. — Comparison of pre-explosion (open circles) and post- 
explosion (filled triangles) cluster SED data. Different lines show 
the change of the SEE) if the SN explosion destroyed a yellow su- 
pergiant (YSG) or a red supergiant (RSG) with a mass indicated 
by the labels. 

the cluster age of ~8-10 Myr. This is consistent with 
the age estimates of S96 found in the previous sections. 

4. DISCUSSION 

In §3.1, §3.2, and §3.3, constraints on the age of S96 
were derived with different techniques. The fitting of 
theoretical SEDs (§3.1) gave possible ages distributed 
between 8 Myr and 40 Myr depending on the cluster 
metallicity and the models applied. The most probable 
solutions turned out to be between 10 and 25 Myr. 
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The fitting of isochrones to the CMDs of the resolved 
stellar population in the vicinity of S96 (§3.2) resulted 
in two distinct populations with ages of 10-16 Myr and 
30-100 Myr. The younger stars seem to be somewhat 
redder, and they are located closer to the cluster center 
than the members of the older population. 

The absence of an H II region emitting in Ha around 
S96 is consistent with the lower age limit of ~ 10 Myr. 
As the simulations with Starburst99 indicate (§3.3), the 
predicted radii of such a cloud after t > 10 Myr would 
decrease below ~ 5 pc, roughly the projected radius of 
S96. 

How can we explain the existence of populations with 
two different ages within such a compact cluster? The 
most likely hypothesis is the capture of field stars by the 
massive ste llar cluster during its formation, as di scussed 
recently bv lPflamrn- Altenburg fe Kroupal (|2007al ) for ex- 
plaining the existence of stars with t « 10-18 Myr within 
the Orion Nebula cluster, where most stars have t < 3 
Myr. The age discrepancy is similar to the case of S96, 
but otherwise the situation is different, because S96 is 
much more massive than the Orion Nebula cluster, and 
the older population resides in the outer region of S96. 

Following the argu ment of 

iPflamm- Altenburg fe Kroupal (|2007af ). the collaps- 
ing pre-cluster cloud may capture nearby field stars 
during its collapse time, which is roughly equal to its 
free-fall timescale, r// « (R^/GMc) 1 ^ 2 , where R c is the 
initial radius of the cloud at the start of the collapse and 
M c is the total mass of the cloud. Adopting R c w 15-20 
pc and M c w (25-100) xlO 3 M Q , the estimated collapse 
time of S96 is r c w 2-8 Myr. Note that this is an 
upper limit for the collapse time, since the mass of the 
cloud can be significantly higher than the final stellar 
mass of the cluster (which was used to estimate M c ). 
Assuming the number densit y of sta rs as 2.44 pc , 
IPflamm- Altenburg fc Kroupal (|2007al ) calculated the 
number of captured field stars within 2.5 pc of the Orion 
Nebula cluster center as 100 < N < 1000 if r c > 2 Myr. 
These results suggest that for S96, which is an order of 
magnitude more massive than the Orion Nebula cluster, 
the number of captured field stars should be substantial, 
even if the number density of the surrounding stars 
is lower than that around the Orion Nebula cluster. 
Moreover, we have studied the stellar content within 
~ 15 pc from the center of S96, instead of 2.5 pc, which 
may also increase the number of captured stars. 

If we interpret the higher age of the stars in the outer 
part of S96 as a result of field-star capture, the pho- 
tometric data suggest that most of the stars resolved 
by ACS are captured, and the young (t ss 10-16 Myr), 
brightest, most massive stars that mostly determine the 
shape of the integrated SED reside within the inner 
(JJ«3 pc), unresolved cluster core. This configuration 
is roughly consistent with that of other young clusters, as 
the brightest, most massive members are generally found 
closest to the center. 

There might be other mechanisms responsible for the 
age dispersion within clusters. A possible hypothesis 
could be continuous star formation within a T c k, 60 Myr 
cloud. This scenario can certainly be ruled out, because 
the SED of such a stellar population is not compatible 
with the observations (see Fig. [TT|) . However, as also sug- 



gested by the referee, there might be a rapid initial star- 
burst that occurred within a region having continuous 
star formation. In this case the resulting SED would be 
a flux-weighted combination of the SEDs from a contin- 
uous SFR (Fig.rjT]) and from a starburst (Fig. [8]). While 
this scenario would certainly be worth study in detail, 
the construction of such customized models is beyond 
t he scope of this paper. 

Pfla mm- Altenburg fc Kroupal (|2008l ) discuss yet an- 
other possibility: gas accretion (and subsequent star for- 
mation) from the nearby interstellar medium (ISM) by a 
massive cluster. However, according to their simulations, 
this process is expected to work only for M c > 10 6 M© 
cluster masses and have a characteristic timescale of a 
few Gyr. Thus, it is probably insignificant for S96, for 
which both the cluster mass and the considered timescale 
are an order of magnitude less. 

The age of S96 is a key parameter in constraining the 
mass of the progenitor of SN 2004dj. The classical the- 
oretical lower limit for the collapse of a stellar core is 
~ 8 Mq, but this can be 1 - 2 Mq smaller depend- 
ing on the treatmen t of core convective overshooting 
(jWooslev et al.l 120021 ). Recent direct identifications of 
Type II-P supe rnova progenitors typically have m asses 
of - 8-15 M Q dMaund et all 120051: iLi et al.l [2006h . and 
none of them clearly exceed M « 20 M Q . 

The fact that SN 2004dj occurred close to the projected 
center of S96 (§2.1.3 and Fig.S]) strongly suggests that its 
progenitor was indeed a cluster member. Although S96 
may contain a significant number of older stars captured 
from the field, it is more probable that a M > 7 Mq star 
is formed during or after the collapse of the pre-cluster 
cloud. Assuming this scenario, the M > 7 Mq limit 
implies t < 60 Myr as an upper limit for the cluster, ac- 
cording to Padova isochrones. This is in good agreement 
with the ages of most of the resolved cluster stars inferred 
from isochrones, because even the members of the older 
population have ages comparable to or less than 60 Myr. 

In Figure[l9]the masses of M > 7 Mq stars are plotted 
as a function of age from the same Padova evolutionary 
tracks as above. The final ages of the curves correspond 
to the last theoretical model for a given initial mass. Note 
that the Padova evolutionary tracks do not extend up to 
the actual moment of core collapse, so the final ages for 
all masses are only lower limits, but an age excess as 
large as ~ 10% is hardly expected. If we accept the ~ 10 
Myr age for S96 as a lower limit inferred from both SED 
fitting and isochrones, this would imply Af prog w 20 Mq 
for the initial mass of the progenitor. 

From fitting the pre -explosion SEP , 

IMaiz-Apellaniz et all (|2004h and IWang et all (|2005l ) 
estimated M prog ~ 12-15 Mq, which agrees very well 
with the most probable age of 1 0-20 Myr fo u nd in 
§3.1 and 3.2. On the other hand, IVinko et all (|2006l ) 
obtained a significantly lower age and higher progenitor 
mass (> 20 Mq) from nearl y the same observed data as 
IMaiz-Apellaniz et al.l (|2004ft . but using different model 
SEDs. This would require T c « 8 Myr, which is lower 
than most of the age estimates discussed above, but may 
not be ruled out entirely, because certain SED models 
indeed predict such young age. However, these earlier 
results were more affected by the age-reddening degen- 
eracy (see §3.1), because of the restricted wavelength 
range of the observed SED. 
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There is yet another way to test the possible mass 
of the progenitor star via the effect of the SN explo- 
sion on the integrated clust er SED, as first suggested by 
iMaiz-Apellaniz et all (|2004D . The explosion of SN 2004dj 
must have changed slightly the supergiant population of 
S96, because one bright (perhaps the brightest) star was 
missing after the SN faded away. This should be appar- 
ent in the cluster SED as well, altering both the over- 
all flux level as well as the spectral shape of the post- 
explosion SED. The difference between the pre- and post- 
explosion SED is approximately the flux spectrum of the 
progenitor star just before explosion. If the progenitor is 
a red supergiant (RSG), then mostly the NIR region of 
the cluster SED will be depressed, while if it is a yellow 
supergiant (YSG), the change will be more pronounced 
in the optical. 

Figure shows a comparison of the observed pre- 
and post-explosion cluster SEDs with the predictions of 
this hypothesis. The lines represent the theoretical post- 
explosion cluster SED if a RSG or a YGS with a given 
mass is removed from the pre-explosion SED. M prog = 
25, 20, and 15 M Q were selected, and their fluxes at 
different bands were determined from the Padova evolu- 
tionary tracks assuming Z = 0.02. 

It is apparent that the removal of a M prog = 25 Mq 
star would cause a strong flux decrease above 5000 A that 
clearly exceeds the uncertainty of the observed flux lev- 
els. Such a massive progenitor is therefore unlikely. The 
lack of a 20 Mq RSG would also result in a similar flux 
depression in the NIR. Unfortunately, no post-explosion 
NIR photometry is at our disposal, so we could not ver- 
ify this prediction. All of the other proposed progenitors 
do not cause the flux to drop significantly below the ob- 
served SED. 

However, as our Keck spectrum suggests (§2.1.2), the 
nebular emission from the SN 2004dj ejecta may also 
have a non-negligible contribution to the observed fluxes 
in the R and / bands via emission lines from Ha, 
[O I] AA6300, 6364, [Fe III A7 155 and possibly [Ca II] 
AA7291,7324 ()Sahu et al.ll2006l ). Also, the lack of post- 
explosion observations in the NIR SED makes the com- 
parison between observations and model predictions un- 
certain at present. More observations, especially in the 
JHK bands, would be very useful to clarify this issue. 
From Figure [20l it seems that the progenitor mass was 
probably < 25 Mq, but a 15 Mq or even a 20 M Q star 
is a possible candidate. 

Putting together all available information, we conclude 
that the new multi-wavelength observations favor a pro- 
genitor star with 12 < Af prog < 20 M Q , if the progenitor 
was a member of the younger population within S96. 
This is consistent with, and perhaps somewhat higher 
than, typical measured SN II-P progenitor masses. How- 
ever, because of the presence of stars belonging to an 
older (~ 60 Myr) population within S96, it cannot be 
ruled out that the progenitor was one of them, which 
would imply M prog rs 7-8 Mq, close to the lower limit 
for such SNe. 

5. CONCLUSIONS 

We have presented late-time photometry of SN 2004dj 
and the surrounding cluster S96, extending the time cov- 
erage of the observational sample up to ~ 1000 d af- 
ter explosion. In the optical, the continuum flux from 



SN 2004dj faded below the level of the integrated flux of 
S96 in 2006 Sep., ~ 800 d after explosion. The pre- and 
post-explosion SEDs of S96 show no significant differ- 
ences in the range 2000-9000 A. The nebular spectrum 
of SN 2004dj at ~ 900 d after explosion was dominated 
by the blue continuum from S96 shortward of 6000 A, 
and by strong Ha, [O I] AA6300, 6363, and [Fe II] A7155 
emission line, characteristic of a typical nebular spectrum 
of a SN II-P. 

We have examined the multi-wavelength observations 
of S96 by different methods, in order to derive constraints 
on the cluster age and evolutionary status. The fitting 
of the cluster SED (using the average of pre- and post- 
explosion fluxes) results in cluster ages distributed be- 
tween ~ 8 and ~ 40 Myr, with the best-fitting solu- 
tions being within 10-20 Myr. The observed reddening 
is E(B — V) « 0.10 ± 0.05 mag; its uncertainty is greatly 
reduced compared with previous studies, due to the in- 
clusion of the UV fluxes from Swift and XMM-Newton. 

S96 appears to be partly resolved in images obtained 
with HST/ ACS on 2005 August 28 (~ 425 d after ex- 
plosion), although the light from SN 2004dj was still 
very strong at that time. We have computed photom- 
etry of the ACS images obtained through the F435IU, 
F606W, and F81AW filters, and combined the magni- 
tudes of the detected stellar sources in color-magnitude 
diagrams. Theoretical isochrones fitted to the observed 
CMDs reveal that the resolved stars in the outskirts of 
the cluster have a bimodal age distribution. The younger 
population consists of stars with ages of 10 < t < 16 Myr, 
while the members of the older one have 30 < t < 100 
Myr. The ages of the older population has a distribution 
that is similar to that of the field stars, not associated 
with S96. This similarity may suggest that about half 
of the cluster stars resolved by ACS were captured from 
the field population during the formation of S96. 

The absence of a visible Ha-emitting cloud around S96 
implies a lower limit for the cluster age of ~8-10 Myr, 
in agreement with the other age estimates. 

The 10 Myr age of S96 would imply a SN 2004dj pro- 
genitor mass of Mp r0 g ~ 20 Mq, while the mass limit for 
core collapse (7-8 Mq) would mean i w 60 Myr for the 
age of the progenitor. This latter limit is consistent with 
the age of the older population within S96, leaving the 
possibility of a low-mass progenitor open. The age of the 
younger population (10-16 Myr) corresponds to M prog « 
12-15 M , which seems to be the most probable mass es- 
timate at present. We verified that even a 20 Mq progen- 
itor would be consistent with the unobservable flux dif- 
ference between the pre- and post-explosion SEDs. How- 
ever, more observations, especially in the JHK bands, 
would be essential to narrow the mass range of the pro- 
genitor. 
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